Ring dark solitons and vortex necklaces in Bose-Einstein condensates 



m 

O 

o 

(N 



O 



I 

o 
o 



> 

O 

(N 
O 
m 
o 



o 



G. Theocharis^, D.J. Frantzeskakis^, P.G. Kevrekidis^, B.A. Malomed'^ and Yuri S. Kivshar* 
^ Department of Physics, University of Athens, Panepistimiopolis, Zografos, Athens 15784, Greece 
^ Department of Mathematics and Statistics, University of Massachusetts, Amherst MA 01003-4515, USA 
Department of Interdisciplinary Studies, Faculty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel 
^ Nonlinear Physics Group, Research School of Physical Sciences and Engineering, 
Australian National University, Canberra AGT 0200, Australia 

We introduce the concept of ring dark solitons in Bose-Einstein condensates. We show that 
relatively shallow rings are not subject to the snake instability, but a deeper ring splits into a robust 
ring-like cluster of vortex pairs, which performs oscillations in the radial and azimuthal directions, 
following the dynamics of the original ring soliton. 



Intensive studies of Bosc-Einstcin condensates (BECs) 
have drawn much attention to the dynamics of non- 
linear excitations such as bright 0] and dark H, ^ |^ soli- 
tons. In particular, dark soUtons in BECs were studied 
in detail and it was found that they are subject to dy- 
namical and thermal instabilities . The experimentally 
observed dynamical instability is due to their quasi- 
ID character: when embedded in a higher dimension, a 
dark-soliton stripe becomes unstable against transverse 
snaking ^,^J. 

Transverse perturbations of dark solitons were exten- 
sively studied in nonlinear optics (where such solutions 
exist in the case of self-defocusing nonlinearity) . Since 
the wavenumber band of the snake instability is limited 
by a maximum wavenumber Qm: the instability can be 
suppressed by bending a soliton stripe to close it into 
an annulus of length L < 27r/Q„i. The resulting ring 
dark solitons (RDS's, i.e., annular troughs on a uniform 
background), first introduced in Ref. 10], were stud- 
ied in optics theoretically 11] and experimentally [T^ . 
We note in passing that bright solitons are unstable to 
collapse in higher dimensions 9], which also pertains to 
bright ring-shaped structures. 

In this Letter we introduce the concept of RDS in 
BEG, as a novel class of solitons which can be experi- 
mentally created by means of known phase-engineering 
techniques '121 ll^ . A principal difference from optics is 
that the RDS dynamics in BEG is temporal, while in op- 
tical media it is of the spatial type [l^- Physical means 
available to control BEG, such as dc and ac magnetic 
fields, are also completely different from those employed 
in optics. Using the perturbation theory |T3 ] and simu- 
lations, we demonstrate that shallow RDS's in BEG are 
long-lived objects, that may be observed experimentally 
on a relevant time scale. On the contrary, deep RDS's 
are subject to snake instability, splitting into ring-shaped 
vortex arrays ( "vortex necklaces" ) which, eventually, re- 
duce to four vortex-antivortex pairs, which perform ro- 
bust double-oscillatory motion in radial and azimuthal 
directions. All these dynamical features are drastically 
different from those known in optics, showing that the 
concept of RDS comprises a much broader range of be- 
haviors than it was known previously. Moreover, we es- 
tablish a link of RDS in BEG to ring fluxons in large-area 



Josephson junctions flUITol. 

The evolution of the BEG is governed by the Gross- 
Pitaevskii equation with a trapping potential V{r) 0]. 
We consider a disk-shaped trap of the form V{r, z) = 
m(w^r^ -I- w^z^)/2, where r"^ = x'^ + y'^ , m is the atom 
mass, LOr^z are the confinement frequencies in the radial 
and axial directions, and f2 = Wr/ujz 1- Then, follow- 
ing Refs. [T^ . one can derive an equation for a normal- 
ized mean-field wave function u{t, r): 



iUt = -(1/2) {Urr + r ^Ur) + \u\'^U + (l/2)f7 



(1) 



We now seek for solutions to Eq. ^ describing rings 
of lower density on a background, which is described by 
the Thomas- Fermi (TF) approximation. For Eq. 
the latter is uq = \/ \i — (1/2) fi^r^ exp {—i^it), where [i 
is the chemical potential. As 17 is small, we can define 
a region where the trapping potential is much smaller 
than \i, then UQ{r,t) « (1/4\/m)(^^)^] exp(— z/xt). 

To describe the dynamics of RDS on top of the back- 
ground Wo, we look for a solution of Eq. QJ in the form 
u = uo{r, t)v{r, t), where the complex field v{r, t) will in- 
troduce the ring soliton. For ^ 1, the most interesting 
case is when the radius of the ring is large enough, so 
that 1/r = 0(ri). In this case, upon redefining t fit, 
r — > ^JJLr, Eq. leads to an effective perturbed nonlin- 
ear Schrodinger (NLS) equation, 



iVt + {l/2)Vrr-{\v\^ ~1)V = P{V), 

where P stands for the effective perturbation, 



(2) 



P{v) = /i-i [(1 - \v\^) vW{r) + {l/2)W'{r)vr - ^{2r)~^Vr\ 

with W{r) = (rir)^/2, all terms in the perturbation P 
being on the same order of smallness. 

We apply the perturbation theory for dark solitons [l3| 
to Eq. 0. We start with the unperturbed dark soliton 
and seek for a ring-like solution to Eq. (01 as v{r,t) = 
cos ^p{t) ■ tanh^ + i sin (p{t), where ^ = cos if (t) [r — R{t)], 
and (p{t) and R{t) are slowly varying phase {\(p\ < 7r/2) 
and radius of the ring soliton. It is straightforward to 
derive perturbation-induced evolution equations: 



dip cos ip dW 



cos Lp 



2n dR 3y/JIR' 



dR 
'dt 



sm ip . 



(3) 
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Combining these, we arrive at an equation of motion for 
the RDS radius: 



1 dW{R) 

2 dR 



1 

3^ 



(4) 



in which we set = 1, as /i can be ehminated from Eqs. 
Q by the transformation t ^/Jlt, ft y/J^^- 

In the hmiting case of a plane soliton, R oo, and 
with cos (fi « 1, it is readily observed that Eq. re- 
duces to an equation of motion for the soliton's radius, 
d^R/dt^ + (il/2)i? = 0, which recovers a known result 
for a quasi-lD dark soliton in a parabolic potential ,6]: it 
oscillates in an harmonic trap with the frequency rJ/\/2- 
On the other hand, in the absence of the trapping poten- 
tial V and for an almost black (deepest) soliton, Eq. Q 
demonstrates that the curvature-induced effective poten- 
tial is U — —(1/3) In i?, which recovers a result known in 
the context of nonlinear optics ^3 • In the present case, a 
combination of the trapping potential and ring curvature 
gives rise to an effective potential well for the soliton's ra- 
dial degree of freedom, U{R) = {l/2){flR)^ - (1/3) Ini?, 
which resembles oscillations of a circular sine- Gordon 
(sG) kink in an axially symmetric anti-trap potential, 
which is possible in large-area Josephson junctions 
(without the anti-trap potential, the circular sG kink 
periodically collapses and bounces back, forming an ex- 
tremely robust pulsating object J^j) 

The above consideration shows that RDS's can be 
found in BEG both as oscillating rings and stationary 
ones, trapped at the bottom of the potential well n(i?), 
with the radius Rq = Q^^y^2/3. For the oscilla- 



i.e 

tory states, the points -Rmin and R between which 
R{t) oscillates can be found, using Eqs. ^) to ehminate 

■ 2 

sm ip: 



R„ 



i2/3)wiO,7j)f^n-\R,, 



where w{k,r]) is the Lambert's w-function defined 
as the inverse of tjiw) = weyi'p{w) 18], the inte- 
ger k is the branch number of the function (k ~ 
corresponds to the principal value), and rj = 
-3W[R{0)] cos^ <f{0) exp {-3W[R{0)]}. 

The possible existence of the stationary and oscillat- 
ing ring solitons is specific to BEGs, where they are sup- 
ported by the trapping potential, while their counterparts 
in nonlinear optics expand indefinitely [Tol |. Stability of 
the ring solitons, trapped at or around R = Rq, against 
transverse perturbations should be tested in direct sim- 
ulations. 

Using the split-step Fourier method ; integrated 
Eq. numerically, with an initial configuration (IG) of 
the form 

u{r, 0) = (1 - n^r^/A) [cos ^(0) tanh Z{r) + sin ip{0)] , 

where Z{r) = (r - Ro)cos(p{0), Q = 0.028, Rq = 28.9, 
and cos (p{0) is the depth of the input soliton. The cases 



FIG. 1: Evolution of a ring dark soliton (RDS) with R{0) = 
28.9 and cos ^(0) = 0.6. (a,b) The initial profile shown by 
cross-section and gray-scale density plots, where RDS corre- 
sponds to a grey ring, (c) RDS shrinks to the minimum radius 
(at t — 40) . (d) Beginning of the emission of dark concentric 
rings (at t — 160). 



of oscillating and stationary RDS can be considered, tak- 
ing cos(^(0) 7^ 1 and cos(^(0) — 1, respectively. Simula- 
tions verify that both oscillating and stationary RDS's 
exist, and their dynamics can be effectively described 
by Eq. Q), up to a certain time. Then, instabilities 
develop: RDS either slowly decays into radiation [for 
cos(/3(0) < 0.67], or, for cos(^(0) > 0.67, snaking sets in, 
leading to formation of vortex-antivortex pairs arranged 
in a robust ring-shaped array (vortex cluster). 

To illustrate these generic scenarios, we first take the 
case with cos(/3(0) = 0.6 and sin(^(0) = —0.8. The cor- 
responding initial structure is shown in Figs. ^a,b). 
According to the analytical results, in this case RDS is 
expected to oscillate with the period T — 240 between 
widely different limits, i?min = 3.8 and i?niax = 69.7, the 
latter being almost at the rim of the BEG cloud, whose 
TF radius isw 70. It is indeed observed that RDS ini- 
tially shrinks, attains the maximum contrast at i? = i?min 
[Fig. n(c)], and bounces back. After reaching i?inax and 
bouncing from it, RDS starts to emit radiation in the 
form of shallow concentric dark rings, as shown in Fig. 
n^d). Due to the radiation loss, RDS becomes shallower 
and, as a result, it accelerates, decreasing the period of 
the oscillations. We observe that RDS performs at least 
three complete cycles of the oscillations before final de- 
cay, which occurs at t w 400. Qualitatively, this dynami- 
cal instability resembles that of a stripe (rectilinear) dark 
soliton in BEGs 2]- 

j^-To translate the results into units relevant to the ex- 
periment 0, we assume ®^Rb condensate of radius 
30/im, containing 20, 000 atoms in a disk-shaped trap 
with ujr = 2Trx 18 Hz and = 27r x 628 Hz. In this case, 
the RDS considered above has the radius Rq = 12.4/im, 
it starts to emit radiation at i ~ 40 ms, and finally de- 
cays at t = 100 ms. This time scale is much larger than 
the hfetime of the dark stripe observed in Refs. 0, |1], 
hence moderately shallow RDS's can be observed too. 

Deep RDS's develop the snake instability, which re- 
sults in the formation of vortex pairs in multiples of 
four, namely 4 [for 0.67 < cosv3(0) < 0.8], 8 [for 
0.8 < cos(^(0) < 0.9], 12 [for 0.9 < costp{0) < 0.95], 
or 16 [for 0.95 < cos(p(0) < 1]. Originally, all the pairs 
are set along a single ring, creating a necklace-like struc- 
ture. The subsequent evolution of the necklaces is char- 
acterized by a transient stage, when quartets of pairs are 
successively expelled off the necklace, drift inward to the 
center of the condensate and disappear there. Eventu- 
ally, there remains a pattern consisting of precisely four 
vortex pairs. They are arranged along a ring that slowly 



oscillates between i?niin and R^, 



i.e., the same lim- 



3 



FIG. 2: Evolution of RDS with R{0) = 28.9 and cos ip{0) = 
0.76. (a) Spontaneous undulations developed on the ring at 
t — 60. (b) Four vortex pairs are formed, at t = 70, as a 
result of the snake instability. 

FIG. 3: Evolution of four vortex pairs created by the insta- 
bility of RDS from Fig. Hfor: (a) t = 120, (b) t = 240, (c) 
t = 420, and (d) t = 540. The vortex pairs remain on the 
ring and move along it, so that the configuration oscillates 
between the x- and +- like configurations, while the ring itself 
periodically shrinks and expands between i?min and R max. 

its between which the initial RDS oscillated prior to the 
onset of the instability. Simultaneously, the vortices and 
antivorticcs perform an oscillatory motion along the ring, 
so that the configuration periodically switches between x- 
and +- like shapes. 

The robust necklace patterns consisting of vortex pairs 
resemble stable clusters of globally linked vortices (of one 
sign, rather than of the vortex-antivortex type) that were 
recently found in a 2D BEC model [23|; however, the 
number of vortices in those clusters could be arbitrary (at 
least, 2, 4, and 8). Another similar object are necklace 
soliton clusters in nonlinear optics, which, however, are 
not stationary, gradually expanding |2l| or rotating . 

The double-oscillatory state persists for long times, 
typically up to i ~ 2000 (which is ~ 500 ms for the 
typical case specified above); still later, due to signifi- 
cant distortion of the condensate as a whole, the radial 
symmetry of the system breaks up, resulting in eventual 
annihilation of all the vortex pairs. To illustrate these 
scenarios, we display two cases, which correspond to situ- 
ations where the instability initially creates the minimum 
(4) or maximum (16) number of vortex pairs. 

First, we consider RDS with cosip{0) = 0.76 and 
^min = 8, -Rmax ~ 58. This IC is vcry similar to that 
in Fig. n^b). It initially shrinks and attains R = i?niin 
[as in Fig. 1(c)]. After bouncing and subsequently ex- 
panding to the rim of the BEC cloud, it starts snaking, 
see Fig. |2Ia), which is a precursor of splitting. Finally, it 
splits into four vortex pairs, see Fig. 2(b). The persistent 
quartet of the pairs arranges itself in a ring configuration. 
The ring performs slow radial oscillations between i?min 
and i?max with a period T 400, which corresponds to 
~ 100 ms, up to t ~ 2000. Simultaneously, the vortices 
and antivorticcs move along the ring, so that they form 
an x-like configuration at i? = i?max [Fig- IS1(£^)]; then 
an octagon [Fig. Ofb)], and then a -I— like pattern. As is 
seen in Fig. |SIc), the latter one shrinks to i? = -Rmin, 
then it bounces and expands, attaining R = i?max [Fig. 
in^d)], and evolves into the x-like pattern, and then the 
cycle repeats itself. 

Finally, we consider the evolution of a black ring soli- 
ton, with cos(/3(0) = 1, which, according to the analytical 
prediction, is expected to be stationary [the correspond- 
ing initial state looks similar to that shown in Fig. ^b)]. 
First, this configuration indeed remains stationary. How- 



FIG. 4: Evolution of an original stationary ring soliton with 
cos ((3(0) — 1 and i?(0) = Ro = 28.9. The initial configuration, 
which is very similar to that in Fig. 0^b), remains undistorted 
up to t = 40. (a) Undulations of RDS at t = 40. (b) Sixteen 
vortex pairs, formed by t = 60 as a result of the snake insta- 
bility. Further snapshots are shown for (c) t = 240 and (d) 
t = 760. The initial set of 16 vortex pairs reduces first to 12, 
then to 8, and eventually to 4 pairs, through successive anni- 
hilations of two quartets at the center of the condensate. Four 
surviving vortex pairs stay at r = _Ro, oscillating between the 
X- and +-like patterns. 



ever. Fig. EJa) shows that, at i = 40 (~ 10 ms), the ring 
starts to snake, which ends up with formation of a neck- 
lace array of sixteen vortex pairs along the ring R = Rq, 
see Fig. ^b). The subsequent evolution of the necklace 
results in annihilation of eight pairs, which occurs in two 
steps. At first, four pairs drift inward, where they dis- 
appear [Fig. EIc)], leaving a nearly rectangular array of 
twelve vortex pairs. Next, the 12-pair pattern expunges 
two quartets of vortex pairs. One quartet again moves in- 
ward and disappears near the center, the other one drifts 
outward, while four vortex pairs stay at R = Rq, see 
Fig. Efc). Then, the four outward- moving vortex pairs 
bounce from the rim of the condensate and move back 
inward, past the quartet that stays put at R — Rq, and 
eventually disappear at the center. Thus, there remains a 
pattern consisting of four vortex pairs, which still reside 
at i? « Rq, see Fig. 4(d) (due to the overall distortion of 
the condensate at this late stage of the evolution, Rq was 
properly adjusted). The vortices and antivorticcs from 
these pairs move along the ring, so that the configura- 
tion performs very slow oscillations between the x- and 
+-like shapes. We have observed almost three complete 
cycles of such oscillations with the period T « 500 (~ 
125 ms). 

In conclusion, we have introduced the concept of ring 
dark solitons in Bose-Einstein condensates, and predicted 
the existence of both oscillatory and stationary solitons. 
Simulations show that perturbation theory accurately de- 
scribes the unperturbed RDS dynamics. However, insta- 
bilities gradually set in and, as a result, shallow RDS's 
slowly decay, while deeper ones develop the snake insta- 
bility. In the latter necklace array consisting of 
vortex-antivortex pairs appears, the number of pairs be- 
ing a multiple of 4. Eventually, it relaxes to a set of four 
pairs which sit on a ring oscillating in the radial direction 
between the same limits which confined the oscillations 
of the original RDS; simultaneously, the pairs perform 
oscillatory motion along the ring. 
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